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1 Introduction 



It is important to develop a string theoretic framework for studying cosmology. The ultimate 
goal of this task is to determine whether string theory can describe basic features of our 
Universe. Despite considerable effort towards this direction over the last few years (see for 
example [1] - [10]), still very little is known about the dynamics of string theory in time- 
dependent, cosmological settings. The purpose of this work is to provide a new class of 
non-trivial string theory cosmological solutions, where some of the difficult issues can be 
explored and analyzed concretely. 

At the classical string level, it seems difficult to obtain exact cosmological solutions [8]. 
Indeed, after extensive studies in the framework of superstring compactifications (with or 
without fluxes), the obtained results appear to be unsuitable for cosmology. In most cases, 
the classical ground states correspond to static Anti-de Sitter or flat backgrounds but not to 
cosmological ones. The same situation appears to be true in the effective supergravity theo- 
ries. Naively, the results obtained in this direction lead to the conclusion that cosmological 
ground states are unlikely to be found in superstring theory. 

Prom our viewpoint this conclusion cannot be correct for two reasons: 

• The first follows from the fact that already exact (to all orders in a') cosmological 
solutions exist, which are described by a two dimensional worldsheet conformal field theory 
based on a gauged Wess-Zumino-Witten model at negative level —\k\: — x M, [2-4]. 

• The second is that quantum and thermal corrections are neglected in the classical 
string/ supergravity regime. 

The first class of stringy cosmological models was studied recently in [5] , where it was shown 
how to define a normalizable wave-function for this class of backgrounds, realizing the Hartle- 
Hawking no-boundary proposal [11] in string theory. Explicit calculable examples were given 
for small values of the level |A;|. As it was shown in [5], these models are intrinsically thermal 
with a temperature below but still close to the Hagedorn temperature. The disadvantage 
of small level \k\, however, is the absence of a semi-classical limit with \k\ arbitrarily large, 
which prevents us from obtaining a clean geometrical picture and studying issues such as 
back-reaction and particle production in a straightforward way. 

Another direction consists of studying the "quantum and thermal cosmological solutions," 
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which are generated dynamically at the quantum level of string theory [6, 12]. Although this 
study looks to be hopeless and out of any systematic control, it turns out that in certain 
cases the quantum and thermal corrections are under control thanks to the special structure 
of the underlying effective supergravity theory in its spontaneously broken supersymmetric 
phase. An effective field theory study has already been initiated in [6, 12]. (See also [13, 14].) 

In order to see how cosmological solutions arise naturally in this context, consider the case 
of a supersymmetric flat string background. At finite temperature the thermal fluctuations 
produce a non-zero energy density that is calculable perturbatively at the full string level. 
The back-reaction on the space-time metric and on certain of the moduli fields gives rise to 
a specific cosmological evolution. For temperatures below the Hagedorn temperature, the 
evolution of the universe is known to be radiation dominated. (See for instance [15, 16] for 
some earlier work in this case and [16] for a review on string gas cosmology.) 

More interesting cases are those where space-time supersymmetry is spontaneously broken 
at the string level either by geometrical [17] or non-geometrical fluxes . In the case where the 
geometrical fluxes are generated via freely acting orbifolds [18] - [23], the stringy quantum 
corrections are under control in a very similar way as the thermal ones. The back-reaction 
of the quantum and thermal corrections on the space-time metric and the moduli fields 
results in deferent kinds of cosmologies depending on the initial amount of supersymmetry 
(Ar = 4,A^ = 2,A^ = 1). 

In this work we restrict attention to four-dimensional backgrounds with initial = 4 or 
N — 2 space-time supersymmetry, obtained by toroidal compactification of the heterotic 
superstring on and T^/Z2-orbifolds. The spontaneous breaking of supersymmetry is im- 
plemented via freely acting orbifolds (as in [18] - [23]). The quantum and thermal corrections 
are determined simultaneously by considering the Euclidean version of the model where all 
coordinates are compactified: x (for the four-dimensional space-time part) xM^ (for 
the internal manifold). Apart from being interesting in their own right, these examples may 
give us useful hinds on how to handle the phenomenologically more relevant N — 1 cases. 
The A" = 1 cases will be studied elsewhere. 

The thermal corrections are implemented by introducing a coupling of the space-time fermion 
number to the string momentum and winding numbers associated to the Euclidean time 
cycle S^. The breaking of supersymmetry is generated by a similar coupling of an internal 
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i?-symmetry charge Qr to the momentum and winding numbers associated to an internal 
spatial cycle Sl^, e.g. the coordinate cycle. 

We stress here, that the thermal and supersymmetry breaking couplings correspond to string 
theoretic generalizations of Scherk-Schwarz compactifications. Two very special mass scales 
appear both associated with the breaking of supersymmetry. These are the temperature 
scale T ~ l/(27ri?o) and the supersymmetry breaking scale M ~ 1/ {27iR<^), with Ro and R5 
the radii of the Euclidean time cycle, S^, and of the internal spatial cycle, S}^, respectively. 
The initially degenerate mass levels of bosons and fermions split by an amount proportional 
to T or M, according to the charges Qf and Qr. This mass sphtting is the signal of su- 
persymmetry breaking and gives rise to a non-trivial free energy density, which incorporates 
simultaneously the thermal corrections and quantum corrections due to the supersymmetry 
breaking boundary conditions along the spatial cycle Sl^. 

At weak coupling, the free energy density can be obtained from the one-loop Euclidean 
string partition function [20] - [22]. The perturbative string amplitudes are free of the usual 
ultraviolet ambiguities that plague a field theoretic approach towards quantum gravity and 
cosmology. For large enough Rq, R^, the Euclidean system is also free of tachyons - the 
presence of tachyons would correspond to infrared instabilities, driving the system towards 
a phase transition [21] - [23]. Therefore, the corresponding energy density and pressure can 
be determined unambiguously, and we can use them as sources in Einstein's equations to 
obtain non-trivial cosmological solutions. This perturbative approach breaks down near the 
initial space-like singularity. We speculate whether this breakdown of perturbation theory 
can be associated with an early universe phase transition. 

The paper is organized as follows. Section 2 is mainly a review, where we also fix most of our 
notations and conventions. We first consider the four-dimensional heterotic string models at 
finite temperature. We obtain the one-loop thermal partition function at the full string level, 
and then we discuss the effective field theory limit at large radius i?o- We also review the 
analogous computation of the one loop string partition function at zero temperature and in 
the case where Susy-breaking boundary conditions arc placed along the internal spatial cycle 
S\j, [18] - [23]. In the large radius limit, the Einstein frame effective potential is proportional 
to the fourth power of the gravitino mass scale, and it can be positive or negative depending 
on the choice of the Susy-breaking operator Qr. 
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In section 3, we consider the case where thermal and quantum corrections due to the su- 
persymmetry breaking are present simultaneously. For the simplest choice Qr — Qf, the 
corresponding one-loop string partition function is invariant under the T M exchange, 
manifesting the underlying temperature/gravitino mass scale duality of the models. This 
duahty is broken by the other allowable choices for the Susy-breaking operator Qji, which 
we classify for both the N — A and the N — 2 orbifold cases. 

In the large radii Rq, limit, the pressure consists of two pieces: the purely thermal part 
which scales as n^T^, with the coefficient being the number of all massless boson/fermion 
pairs in the initially supersymmetric theory, and another potential-like piece which scales 
as UyM'^ and with the coefficient riy being positive or negative depending on the choice of 
the operator Qji. In both pieces, the rest of the dependence on the scales T and M can be 
expressed neatly in terms of non-holomorphic Eisenstein series of order 5/2 whose variable 
is the complex structure-hke ratio M/T. In addition, we incorporate the effects of small, 
continuous Wilson line deformations in our computation. Wilson lines along any of the 
internal spatial cycles, other than Sj^, introduce new mass scales, and pieces proportional 
to ~ and ~ arise in the effective thermodynamic quantities. 

In section 4 we present our ansatz for the induced cosmological solutions. These are homo- 
geneous and isotropic cosmologies for which the Susy-breaking scales T and M as well as the 
inverse of the scale factor 1/a evolve the same way in time, and so the ratio of any two of 
these quantities is constant. The form of this ansatz is dictated by the scaling properties of 
the effective energy density and pressure. The compatibility of the gravitational field equa- 
tions with the equation of motion of the scalar modulus controlling the size of the gravitino 
mass scale fixes the ratio M/T. By solving the compatibility equations numerically, we find 
that in the absence of Wilson lines along S]^, non-trivial four dimensional solutions exist 
when ny is negative and the ratio l^iyl/n^ is small enough. These conditions are satisfied 
by various models we describe explicitly in the paper. When we include Wilson lines along 
SIj, the value of the ratio M/T for some of the solutions can be large or small, and so we 
can have models with a hierarchy for the scales M and T. 

Having solved the compatibility equations, the time-dependence of the system is governed 
solely by the familiar Friedmann-Hubble equation. There is a radiation term, Cr/a^, whose 
coefficient is positive in our examples. An effective curvature term, —k/a?, can be gen- 
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erated by turning on Wilson line deformations. The sign of k can be a priori positive or 
negative, depending on the model. When we turn on the kinetic terms of some of the extra 
flat moduli, we generate an additional term that scales as c^/a^ (with positive). 

In section 5, we solve the Priedmann-Hubble equation for the various possible cases, and we 
elaborate on the properties of the cosmological solutions: 

• When Cr > 0, we have standard hot big bang cosmologies with an intermediate radiation 
dominated era. The late time behavior is governed by the spatial curvature of the models. 

• We also consider a priori possible exotic models characterized by < 0. A big bang 
occurs when Cm > 0. The cosmological evolution always ends with a big crunch when A; > 0. 
The case k < however is more interesting. It involves either a first or second order phase 
transition between the big bang cosmology and a linearly expanding universe. The first case 
corresponds to a tunneling effect involving a gravitational instanton, while the transition is 
smooth in the second case. If the first order transition does not occur, the universe ends in 
a big crunch. 

We finish with our conclusions and directions for future research. 



2 Thermal and quantum corrections in heterotic back- 
grounds 

Our starting point is the class of four dimensional string backgrounds obtained by toroidal 
compactification of the heterotic string on and T^/'Z,2 orbifolds. Initially the amount of 
space-time supersymmetry is A^4 = 4 for the case of compactification on the torus and 
N4 = 2 for the orbifold compactifications, and the four dimensional space-time metric is flat. 
Space-time supersymmetry is then spontaneously broken by introducing Scherk-Schwarz 
boundary conditions on an internal spatial cycle and/or by thermal corrections. Due to the 
supersymmetry breaking, the one-loop string partition function is non- vanishing, giving rise 
to an effective potential. Our aim is to determine the back-reaction on the initially flat 
metric and moduli flelds. 

At the one-loop level, the four dimensional string frame effective action is given by 

j d'x,/^d^ + +•••)- Vstring) , (2.1) 

where is the Ad dilaton held and the ellipses stand for the kinetic terms of other moduli 
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fields (to be specified later). At zero temperature, the effective potential Vstrmg can be 
obtained from the one-loop Euclidean string partition function as follows: 

TT = -Vstring , (2.2) 
1/4 

with V4 the 4d Euclidean volume. The absence of a dilaton factor multiplying the potential 
term in the action is due to the fact that this arises at the one loop level. 

At finite temperature, the one-loop Euclidean partition function determines the free energy 
density and pressure to this order 

TV — —^String = -Pstring. (2.3) 
1/4 

The subscript indicates that these densities are defined with respect to the string frame met- 
ric. The relevant Euclidean amplitude incorporates simultaneously the thermal corrections 
and quantum corrections which arise from the spontaneous breaking of supersymmetry and 
which are present even at zero temperature. 

In order to determine the back-reaction of the (thermal and/or) quantum corrections, it is 
convenient to work in the Einstein frame where there is no mixing between the metric and 
the dilaton kinetic terms. We define as usual the complex field 5", 

5 = e-2'^ + iX, (2.4) 

where x is the axion field. Then after the Einstein rescahng of the metric, the one loop 
effective action becomes: 



S — J d^x\/ — det g 



- g^" Kjj- d^<^>ld,^j - ^ Vstnng(*/, */) 



(2.5) 



where Kij is the metric on the scalar field manifold {$/}, which is parameterized by various 
compactification moduli including the field S. This manifold includes also the main moduli 
fields T/, Uj, I = 1,2, 3, which are the volume and complex structure moduli of the three 
internal 2-cycles respectively. We notice that in the Einstein frame the effective potential, 
Vstring, is rescaled by a factor where s — Re(S') = e~'^'^. Taking this rescaling into 

account, we have 

VEin = ^"^String. (2.6) 

s 

This relation will be crucial for our work later on. (We will always work in gravitational 
mass units, with Mq = ^g= = 2.4 x 10^^ GeV). 



Keeping only the main moduli fields {S, Tj, Ui}, their kinetic terms are determined in terms 
of the Kalher potential K [24,25]: 

X = -log {S + S)-J2 log iTi + fj)-J2 log (Uj + Uj) (2.7) 

I I 

with Kjj = djdjK. The classical supcrpotential depends on the way supersymmetry is 
broken. Generically string backgrounds with spontaneously broken supersymmetry are fiat 
at the classical level due to the no-scale structure of the effective supergravity theory [25]. 
Once the thermal and/or quantum corrections are taken into account, we obtain in some 
cases interesting cosmological solutions. 

2.1 Heterotic super symmetric backgrounds at finite temperature 

In order to fix our notations and conventions, wc first consider the case of an exact supcrsym- 
metric background at finite temperature [21] - [23]. For definiteness we choose the heterotic 
string with maximal space-time supersymmetry (A^4 — 4). All nine spatial directions as well 
as the Euclidean time are compactified on a ten dimensional torus. At zero temperature, 
the Euclidean string partition function is zero due to space-time supersymmetry. At finite 
temperature however the result is a well defined finite quantity. Indeed, at genus one the 
string partition function is given by: 

7 - / ^'^^^ 1 \a+b+ab n ral4 ^(10.26) [g] 

where r(io,26) [&] is a shifted Narain lattice (which we specify more precisely below). The 
non-vanishing of the partition function is due to the non-trivial coupling of the lattice to 
the spin structures (a, 6). Here, the argument a is zero for space-time bosons and one for 
space-time fermions. The spin/statistics connection and modular invariance require that the 
unshifted r(i^i) sub-lattice of the Euchdean time cycle 

r(,,) ^ ^o(Imr)-^ e-«§^ (2.9) 

m,n 

be replaced as follows: 

r(l,l) ^ ^i?o(ImT)-^ e-'^-^o^'^ gi7r(ma+n6+mn) _ (2.IO) 

m,n 
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Redefining 

m^2m + g, n^2n + h, (2.11) 

where g, h are integers defined modulo 2, and introducing the notation i) [^] for a shifted 
lattice, 

rci. G] = E ^(i--)-^ ^-"^ . (2.12) 

the thermal partition function takes the form: 

f 1 Nffa+fe6+% / Na+6+a6 /J ral4 ^(9.25)r(l,l) [g] 

2 2. (-) (-) U • 

Defining a = a — h and b = b — g and using the Jacobi identity 

(2.14) 

(a,S) 

we obtain 

^-il^i:-(-r"''[i:rr,M.[;];;^^. (2.15, 

The temperature in string frame is given by Tstring = l/(27ri?o). 

Since our aim is the study of induced cosmological solutions in 3 + 1 dimensions, we consider 
the case for which the radii of three spatial directions are very large: — Ry — Rz = R ^ 
In this case the three dimensional spatial volume factorizes 

r(3,3) = R' (Imr)-i = ^ (Imr)-i. (2.16) 
Using the expression for the r(i^i) [^'] shifted lattice we obtain: 

Z — — (27ri?o)V'3 -^String = V4 Pstring 

- _ i ^^^^ V ^.i^2 i2^+.K2„+.).p 4 r(6,22) . . 

(27r)4/; 4ImT3 ^ ^ L1+9J (^)i2 ^(^)24 ' y^-^') 

{n,rn),{h,g) 

where V4 = {2t:Rq)V^ is the four dimensional space- time volume, ^string the free energy 
density and Pstring the pressure in string frame. 

Before we proceed further, we make some comments: 

• The sector (/i, g) = (0, 0) gives zero contribution. This is due to the fact that we started 
with a supersymmetric background. 
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E(-) 



a+b+ab 



d+h 
b+g 



• In the odd winding sector, h = 1, the partition function diverges when Rq is between 
the Hagedorn radius Rh — (a/2 + l)/2 and its dual 1/Rh- < -^o < Rh- The divergence 
is due to a winding state that is tachyonic when Rq takes values in this range, and it signals 
a phase transition around the Hagedorn temperature [21] - [23]. In this paper we study the 
regime Rq > Rh, where there is no tachyon and the odd winding sector is exponentially 
suppressed. The high temperature regime and the cosmological consequences of the phase 
transition will be examined in future work [32]. 

• When Rq ^ 1, the contributions of the oscillator states are also exponentially sup- 
pressed, provided that the moduli parameterizing the internal ^{6,22) lattice are of order 
unity. 

2.2 The effective field theory in the large Rq limit 

As we already mentioned, the h — 1 sector of the theory gives exponentially suppressed con- 
tributions of order O{e~^o). Also, the {h, g) — (0, 0) sector vanishes due to supersymmetry. 
Thus for large i?o, only the sector (/i, g) = (0, 1) contributes significantly. Using the identity: 

r(i,i)(i?o) = r(i,i) [°] + r(i,i) [?] + v^,,,^ [J] + v^,,,-^ [i] (2.18) 

and neglecting the h = 1 sectors, we may replace 

r(i,i) [?] ^ r(i,i)(i?o) - r(i,i) [^] = t^,^,){Rq) - \t(,,,pr^) (2.19) 

in the integral expression for Z . For each lattice term we decompose the contribution in 
modular orbits: {m,n) — (0,0) and {m^n) ^ (0,0). For {m,n) ^ (0,0), the integration over 
the fundamental domain is equivalent with the integration over the whole strip but with 
n — 0. The (0, 0) contribution is integrated over the fundamental domain. Now the (0, 0) 
contribution of r(i^i)(i?o) cancels the one of |r(i^i)(2it!o), and we are left with the integration 
over the whole strip: 

^ f V e-^^o'^ [1]^ ^(6-^^) . (2.20) 

(27r)4 ill 4ImT3 ^ r]{Ty^ r?(f)24 ^ 

The integral over ti imposes the left-right level matching condition. The left-moving part 
contains the ratio 

n rii4 

^ = 2^ + C»(e-'^^2), (2.21) 
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which imphes that the lowest contribution is at the massless level. Thus after the integration 
over Ti {t2 = t), the partition function takes the form 



(2vr; 



--^"'"^ ^0 + E Di,) e-^^) , (2.22) 



where -D(/i) denotes the multiplicity of the mass level fi and 2^ Dq is the multiplicity of the 
massless level. Changing the integration variable by setting t = TiR^{2m + 1)^ x, we have: 

1^4 V- 1 dx 



vr2(27ri?o) 



4 

m 



Now the second term in the parenthesis is exponentially suppressed when the masses /x are 
of order (or close) to the string oscillator mass scale. This will be the case when all of the 
internal radii and the Wilson-line moduli of the r(g 22) lattice are of order unity. For this 
specific case, the partition function simplifies to 

\^{2tiRoY ^ i2m + lY~ 48 (27ri?o)^ " 3 16 'string, l^-^^j 

where n* = 2^Dq is the number of the massless boson/fermion pairs in the theory. The free 
energy density and pressure in string frame are given by 

X Ti TT -^String 

-Pstring = "-^String = " — • (2.25) 



In the Einstein frame, energy densities are rescaled by a factor 1/s^ as in Eq. (2.6). Thus 
the pressure and free energy density in this frame are given by 

„ _ ^ _ 1 ^Itring _ 1 n*TT^ 

-TEin — — ^Ein — 77 77 ^ — 77 777 , (2.26) 

3 16 s"^ 3 16 

where T = T^trmg/ Vs is the proper temperature in the Einstein frame. This result is expected 
from the effective field theory point of view. When only massless states are thermally excited, 
the field theory expression for the pressure is given by 

+ (2.27) 
where ns and np are the numbers of massless bosonic and fermionic degrees of freedom 



respectively. When ub = np = n*, as in a supersymmetric theory, we recover Eq. (2.26) 
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2.3 Spontaneous breaking of supersymmetry at zero temperature 



In this case we consider the same class of heterotic models, but now the breaking of super- 
symmetry arises due to the coupling of the space-time fermion number to the momentum 
and winding quantum numbers of an internal spatial cycle [18] - [23]. Since the temperature 
is taken to be zero, the spin structures (a, b) do not couple to the quantum numbers of the 
Euclidean time cycle which will be taken to be very large. We also consider the case where 
three additional spatial directions are large. Following similar steps to the purely thermal 
case, the partition function is given by 

(n,m),{h,g) 

where now V5 = ¥4(271 R^) is a five dimensional volume and the 21) lattice parameterizes 
the internal space. Here also, the h = 1 sectors give exponentially suppressed contributions 
0{e~^^, and the {h, g) = (0, 0) sector vanishes due to supersymmetry. The rest of the steps 
can be repeated as in the derivation above to find 

K dt ^ „o2(2m+l)2 



"'0 2*2 _ ^ / 



(27r)5 

which after the change of variables t = nR\{2m + 1)^ x gives 

^|2m + l|5 7o 2x1 V ^ ^ / 



(2.29) 



7r2(27ri?5)5 
For /X of order unity, this simplifies to 



(2.30) 



2(1-2-) ^^^rz^-A^ (2.31) 
7r2 



with n* = 2^Dq. 



This result was expected from the effective field theory point of view. Indeed in a theory with 
spontaneously broken N — 4 supersymmetry, the one loop effective potential receives a non- 
zero contribution proportional to the mass super-trace StrM^, which in turn is proportional 
to the fourth power of the gravitino mass. The super-traces StrM."^ vanish lor: n < N — A. 
In the example of supersymmetry breaking we examined above, the masses of the states are 
shifted according to their spin. For initially massless states, the mass after supersymmetry 
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breaking becomes : 

Ml ^. (2.32) 

This shows that the string frame gravitino mass is of order Mstring ~ I/-R5 and thus StrAi^ ~ 
c/R^. Including the contributions from all Kaluza-Klein states, one obtains the result given 



in formula (2.31). We obtain for the string frame effective potential: 



Vst.n, = = -2 (1 - 2-^) n* -4—. (2.33) 



In the Einstein frame, we have VEm = ^Vstring - see Eq. (2.6) - so that 



where ti = Re(Ti), ui = Re([/i), and M = l/^stiUiY^"^ is the gravitino mass scale in the 
Einstein frame. 

We stress here that the one loop effective potential depends only on the gravitino mass scale, 
which in turn depends only on the product of the s, ti and ui moduli. This suggests to freeze 
all moduli and keep only the diagonal combination 

3 log z = log s + log ti + log Ui . (2.35) 

The Kalher potential of the diagonal modulus Z, (with z = Re(Z)), takes the well known 
SU{1, 1) structure [25] 

K = -3\og{Z + Z) . (2.36) 
This gives rise to the kinetic term and gravitino mass scale. 

Freezing \m.Z and defining the field $ by 



one finds the kinetic term 
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The choice o? — 3/2 normahzes canonically the kinetic term of the modulus The potential 
for this particular model is: 

VEin($) = -Cv = -Cv e^"*, « = . (2.40) 

Observe that in this simple model the sign of the potential is negative. As we now explain, 
we can construct models with a positive potential, but with the rest of the dependence on 
the modulus $ being the same. All we have to do is to couple the momentum and winding 
numbers of the Scherk-Schwarz cycle not only to the space-time fermion number but also 
to another internal charge. For example consider the x E'^ heterotic string on and 
instead of coupling just to Qf, we couple to Qp + Qes+ Q'es^ where Qes denotes the charge 
of an Eg representation decomposed in terms of 5*0(16) ones, and similarly for Q'^^- These 
charges take half integer values for the spinorial representations and integer values for the 
others. The initial Susy-breaking co-cycle gets modified as follows 

^_^ag+bh+hg ^ ^_-j(a+7+7')a+(6+5+5')/i+/ip ^2 4]^^ 

where as before the argument a is one for space-time fermions and zero for space-time 
bosons, and (7,7') = (1, 1) for the spinorial representations of SO{16) x 30(16)' and (0,0) 
for the adjoint representations. This operation breaks explicitly the x E^ gauge group to 
5*0(16) X 50(16)'. Proceeding in similar way as in the previous example, one finds: 

^.2(1-2-=) «5£(|)«-pJ|p, (2.42, 

where 

r = 2^ [ [2]x2,3 + [6]r6 + [120 - 128]b3 + [120 - 128]^/ ] = - 2^ x 8 = - 64. (2.43) 

In the previous example only positive signs appear in the above formula since there is no 
coupling of the Scherk-Schwarz lattice quantum numbers to the E^ x E'^ charges, giving the 
value n* = 2^ X 504. The reversing of sign for some representations indicates that it is for the 
bosons that the masses are shifted and not for the fermions in the corresponding multiplet. 

We note that in the N = A case, we cannot change the left-multiplicity since all of the 
left-moving R- charges arc equivalent as required by symmetry. This however is not true 
for the N = 2 and N = 1 cases. Consider for instance the class of = 2 supersymmetric 
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backgrounds obtained by compactifying the heterotic string on a T"^ j%2 orbifold (e.g. the Z2- 
orbifold hmit of the CY-compactification). In this class of models (see for instance [27]) 
four internal supercoordinates are twisted and the corresponding four internal R-charges are 
half-shifted. The Euclidean partition function is given by 



2 



p 4Imr 4 ^ riirY 

r(i,i)(i?o) r(3,3) (space) ^ 

r^(r)2 r7(f )2 [oj ^(4,4+n,) [gJ • (2-44) 

Here Z(2^2+no) is the contribution of two internal coordinate^ (X5, Xg) and ng-right moving 
world-sheet bosons Before supersymmetry breaking, the corresponding (2, 2 + no)-lattice 
is unshifted. Z(4^4+„() stands for the contribution of four internal coordinates (X7, Xg, Xg, Xio) 
all of which are Z2-twisted by (//, G), and nj-right moving world-sheet bosons 0/ which can 
be Z2-twisted breaking part of the initial gauge group. The 6'-function terms come from 
the contribution of the left-moving world-sheet fermions. Four of them are Z2-twisted by 
{H,G). The contribution associated to the space-time bosons is when a = 0, while the one 
associated to the space-time fermions is when a = 1. 

From the above supersymmetric N = 2 partition function, the thermal partition function 
is obtained in a way similar to the X = 4 example, by the following replacement of the 
Euclidean time sub-lattice: 

r(i,i)(i?o) r(i,i) [J;] (Ro) (_)9i«+/^i^+/^i9i. (2.45) 

In the case of Scherk-Schwarz spontaneous supersymmetry breaking, the partition function 
can be obtained by a similar replacement of the internal X5 coordinate lattice, either by 
utilizing the same operator Qp 

r(i,i)(i?5) ^ r(i,i) [J^^] (R,) (-y^-+f^^b+h,g, ^^aq) 

or by utilizing an R-symmetry operator associated to one of the twisted complex planes 

r(i,i)(i?5) r(i,i) (i?5) (-)^'^(-+^)+^^(''+«)+'^^^^ (2.47) 

These are in fact the only two possibilities involving left-moving R charges since all others 
are equivalent choices. However, many other choices exist by utilizing parity-like operators 
^In our notations, the space-time coordinates are ^o....,3i while the internal ones are ^5,....io- 
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involving the right moving gauge charges ^ 7^, as in the exphcit example of 50(16) x 5*0(16)' 
spinorial representations we gave above: 

r(i,i)(i?5) 1(1,1) [5^1 (^5) {-y^^''+"+^^^^+^'^''+^+^~^^^+^'^\ (2.48) 

In the next section we examine representative examples in the case where thermal and 
spontaneous Susy breaking operations are present simultaneously. 



3 Thermal and spontaneous breaking of Susy 

The most interesting situation for cosmological applications is the case where spontaneous 
supersymmetry breaking and thermal corrections are taken into account simultaneously. 



3.1 Untwisted sector 



The untwisted sector of the = 2 case, (if, G) = (0,0) in Eq. ( |2.44| )f1 has an = 4 



structure and thus all choices for the left R-symmetry operators are equivalent. The quantum 
numbers of the Euclidean time cycle and the internal Xs-cycle are coupled to the spin 
structures (a, b) in the same way. After performing the Jacobi theta-function identity the 
partition function becomes: 

(ni,mi),{hi,gi) (n2,m2),(/i2,92) 

^^fl2 |2mi+gi+(2ni+hi)r|2 j^2 \2m2+32 + (2n2+h^).\^ n+hl+hil'^ ^(5,21) lN 

^ mr e '^Ll + <?l+92j ^(-^^12 ^(•^^24 • ^'^■^^ 

The factor of 1/2 is due to the Z2 orbifolding of the A^ = 4 theory. 

Proceeding as in the simpler examples before and neglecting the hi = 1 and h2 = 1 sectors 
for large Rq, R5, the non-zero contributions to the partition function occur when gi+g2 = 1- 
Assuming also that all other moduli are of order unity, the only non-exponentially suppressed 
contributions come from the zero mass left- and right-levels. We obtain 

V, V- (1 - (-)^^+^^) r dt V- _^.2 (2^i+.i) ^ 

^untwist -^0 / 

91,92 mi,m,2 



2'^o E^^^^ X ^ E .-'"i-^-"^-^. (3.2) 



^For hi = h2 = (even windings), the sector {H,G) = (0,1) gives zero net contribution due to the 
identity i E]'^ [?+i] 9 [Li] = 0. 
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which after the change of variables t — n {Rl{2mi + gi) + i?g(2m2 + 52) ) x gives 

_ 4L>or(|) V, ^ 1 

■'untwist — 6 /'o—NR / J 

mi,m2 



(27r)5 ^f- (Rl(2rm + 1)^ + it:i(2m2)^)i 



I r (i) ^ 

(27r)5 {Rl{2rmr + Rl{2m2 + ir)l' 

This expression is symmetric under the Rq <-> R5 exchange. This is suggestive of a temper- 
ature/gravitino mass, T/M, duahty. This duahty will be broken when the supersymmetry 
breaking arises due to the coupling to a different Qr charge than Qp. 

To obtain the effective four dimensional pressure, we must factorize out the space-time 
volume V4. To this extent it is convenient to define the complex structure-like ratio 

Ro M 

= (3.4) 
and re-write the partition function in the following way 



■'untwist 5 



E 



ttI (27ri?o)^ ^ |(2mi + l)ra + 2m2p 

^ ' mi,m2 

I r (i) V, 1 

7ri (27ri?5)^ ^ \2miiu + {2m2 + ^'^ 



mi,m2 

Define the function 



^^^^ = ^ \{2rm + l)iu + 2mr ^'"'^ 
which we can express in terms of Eisenstein functions of order 5/2: 

f{u) = ui {^^E,„ (^1) - j,E,/,{iu)^ , (3.7) 



where 



(m,n)^(0,0) 

Then the pressure in the Einstein frame can be written as 

Puntwist = /(«) + CT ffi/^, (3.9) 

u 

where 

= = <nt (3.10) 
16 



Here n*^^ = ADq is the number of massless boson/fermion pairs in the untwisted sector. It 
is smaller by a factor of 1/2 from the corresponding number in the N — 4 case due to the 
Z2-orbifolding. In this particular model the coefficients C|.°* and Cy^* are equal due to the 
underlying gravitino mass/temperature duality. For fixed u the first term stands for the 
thermal contribution to the pressure while the second term stands for minus the effective 
potential. 

We note that the coefficient C^"* is fixed and positive as it is determined by the number of all 
massless boson/fermion pairs in the untwisted sector of the initially supersymmetric theory: 
= ADq. In general, the coefficient Cy"* will depend on the precise way supersymmetry 
is broken. As we have demonstrated in the previous section, it can take both negative and 
positive values, depending on how the Susy-breaking operator couples to the right movers: 
C™* ~ Thus in general the temperature/gravitino mass scale duahty will be broken. 

Let us discuss the large u limit, which can be obtained by taking to be small (but still 
parametrically larger as compared to the string scale), while taking Rq to be much larger. In 
this limit we expect to find a four dimensional system at finite temperature, for which only 
massless bosonic degrees of freedom are thermally excited. All fermions attain a mass from 
the Scherk-Schwarz boundary conditions along the cycle, and this mass is much bigger 
than the temperature for large u. Therefore they can be integrated out giving a temperature 
independent contribution to the pressure of order l/Rf- Setting u — 1/u, we have in the 
limit u — > (li — > oo): 

f{u) = f{l/u) - E r 77-. T^^r^ = I E (o = ^ (3.11) 



and 



^^„-oo((2mi + l)2 + 4x2). 
f{l/u) 1 1 ^ f°° dx 



u 



= 2(1-2-5) ^(5^+ 1 y 1 



24 ^ mi 

= 2(l-2-V(5) + Y^x^^. (3.12) 



Using these results, we find 

1 2^Do TT^T^ , C(5)r(i) 



untwist 



30 



2 (1 - 2-5) ' s^^^ (4Do) (3.13) 
7r2 
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for the first two leading terms for small u = 1/u. We have used the relation r(|) = Sy/n/A. 
The first contribution arises from the thermally excited massless bosons. As compared to 



Eq. (2.27) with ub = 2^Dq and np = 0, it is off by a factor of 1/2 due to the Z2-orbifolding. 



Similarly, the second term is off by a factor of 1/2 as compared to Eq. (2.31) due to the 



orbifolding. For large u, the potential term is dominant. For generic values of u both 



fermions and bosons contribute to the thermal piece as in Eq. (3.9), with the contribution 
depending on the number of massless states at zero temperature and before the breaking 
of supersymmetry. Finally for small u, the system is essentially a five dimensional purely 
thermal system. 

3.2 Twisted sector, H = 1 

The contributions of the twisted sectors in the large Rq, limit depend on the number of 
the massless twisted states before the supersymmetry breaking, and can be determined in 
a similar way as before. However, there is a class of models where the Z2-orbifolding acts 
freely, without any fixed points, and therefore there are no massless states in the twisted 
sectors. For this class of models, the whole contribution to the one-loop partition function, 
in the large Rq, R^ limit, is that of the massless untwisted sector states we have already 
determined. One example with this property is when the Z2-twists {H, G) are accompanied 
with a shift of the r(i^i)(i?6) sub-lattice. This operation leads to the modification of -2^(2, lo) 



in Eq. (2.44), where we also set {no,nt) = (8,8), as follows: 



7 roi ,7 mi r(i,i)(i?5) 1(1,1) (i^e) [g] 1^ []f , . 



If Rq is sufficiently large, the coupling of the {H, G)-shift of the lattice to the twisted partition 
function .^'(4,12) [q] ensures the absence of massless states in the twisted sector. 

In other situations, there are massless states in the twisted sector. Before the supersymmetry 
breaking, the number of massless bosons is equal to the number of massless fermions with 
a multiplicity nl^^^^. Proceeding as in the untwisted sector, and neglecting the /ii,/i2 = 1 
sectors, one finds that there is only a non zero contribution when Qi + g2 = ^- The relative 
sign of the thermal part ~ and the supersymmetry breaking part ~ M*^ depends on the 
choice of the operators Qp and Qj?. When there is no coupling to the right-moving gauge 
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charges we obtain: 

u 

where 

r (^) 

/^twist * v 2 / 

712 



r C^) 

^twist c ri* V 2 / 

'^V ~ t 'Hwist 5 
7T2 

. ^ ,'_\(Qfl-QF) 



(3.16) 



Here, we have for the coefficient e: 

• e = 1 when Qfj = Qp. 

• e = — 1 when ^ Qf- 



In the later case, (— )('3«"Q-f) = [—)^ = —1, see Eq. (2.47), and H = 1 in the twisted sector. 
The change of sign indicates that it is the bosons that are becoming massive because of the 
supersymmetry breaking. This is related to a mechanism for moduli stabilization induced 
by geometrical fluxes [17]. 

Adding the contributions of the untwisted and twisted sectors together we obtain for the 
pressure 

p= CtT^ f{u) + Cv Zi^^ (3.17) 

u 

where Ct = + C^^^^ and likewise for Cy- The sign of the thermal contribution is always 
positive, 

Ct = (3.18) 

712 

— ^unt + ^twisf The coefficient multiplying the supersymmetry breaking part is given by 

Cy = ^^^l^ffi (3.19) 
712 

with Uy = 72* jj^. + e In general riy can be positive or negative depending on the model. 



3.3 An explicit example 

As an example we consider the Eg x Eg heterotic string on a T^/Z2 orbifold, whose initially 
supersymmetric partition function is obtained by setting 

r(i,i)(P5) r^i,i)m 1 ^ 9 []f 



r]{ry r/(r)2 2 ^ r]{T) 
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and 



"(4,12) 



[S] 



(4,4) 



[g] 



7',<5' 



■y'+H 
S'+G 



S'-G 



(3.21) 



in Eq. (2.44). We shall use an R-symmetry operator associated to one of the twisted complex 
planes for breaking the supersymmetry, replacing the i)(i?5) lattice as in Eq. (2.47). In 
the twisted sectors, {H, G) ^ (0, 0), the internal 4) shifted lattice is given by 



^(4,4) [g] 



1 2 • 



The orbifolding breaks the Es x Es gauge group to EsxEjX SU{2). Under E^ 
the 248-dimensional adjoint representation of E^ decomposes as 



fl,3) 



(56,21 



:i33,i). 



The untwisted sector contains 2'^ x 504 massless states giving the value n* 



unt 



(3.22) 

EtxSU{2), 

(3.23) 
4 X 504 



for the total number of boson/fermion pairs. These numbers arise as follows. In terms of 
world-sheet left / right movers the number of bosonic degrees of freedom is given by 

<„t = 4 X 504 = [4]^2.3,5,6 X ([4]x2,3.5,6 + [248]^, + [133]^,, + [3]5t/(2)) + 

[4]^7,8,9.io X ([4]x7,8,9,io + [2]su(2) X [56]^;,) . (3.24) 

The first line gives the bosonic content of a c? = 6 supergravity multiplet, a tensor multiplet 
and a vector multiplet in the adjoint of the E^ x E^ x SU{2) gauge group. The second line 
gives the bosonic content of four uncharged and one charged hyper-multiplets. The number 
of fermionic degrees of freedom follows by supersymmetry. At finite temperature and when 
Susy is broken, the contribution of the massless untwisted states is determined as before (see 



Eq. (3.5)) 



Next we analyze the contribution of states in the twisted sectors, H = 1. For large RqjR^, 
we may neglect the hi = 1 and /i2 = 1 sectors. Also for Rq of order unity we may set 
r(i^i)(i?6) — 1- Setting H = 1, the partition function becomes in this limit 



^twisted 



drdr 



,2 {2mi+9ir 



2 (2m2+92) 
Imr 



X 



iE( 

{a,b,G) 



(mi,9i),(m2,92) 
a+b+ab ^_-^agi ^_'j(a+l)g2 ^ |^^^^ j ' 
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X 



,2 'E^^Bi' 5E»1i' 



e 



7'+l 
S'+G 



9 



y-i 

5'-G 



. (3.25) 



^ Li+gJ ^ Li+gJ 7,5 7',5' 

The non- vanishing contributions arise when (7i + (72 = 1- Non-exponent ially suppressed 
contributions arise only at the zero mass leveL To obtain them, we expand the integrand in 
powers of g = e^'^*^. 

When (71 + = 1 , we have for the left movers 
1 

(a,6) 

1 



\b+ab 



{a,b) 
32/32/32 

' for G = or 



91 

4(-)^^(l + 0(gi/2)). 



\ 92/32/32 

^ ^2^^ for G=l 



(3.26) 



For the right movers we have 
1 1 



9 



7'+l 
<5'+G 



7'-l 
5'-G 



4j (l + H''4gV^ + ...) (1 + 240 g + ...)[gV^((-)^2 + 56g^/^ + ...)] 



9 



i^[(-)^2 + 8gV2 + 56g^/2 + 0(g)] . 



(3.27) 



When we add the contributions of the {H,G) = (1,0) and {H,G) = (1,1) twisted sectors 
together, we find that the lowest right-mass level is at zero mass as it is the case for the 
lowest left-mass level. 

Using these results we obtain the contribution to the partition function of the massless 
twisted states: 



2^(56 + 8) 



"192 _ (\9i 



(27r)5 ^ 2 

91 >92 



dt 



2t2 



_^ p2 (2"'l+9l)^ ^^r2 (2^2+S29 
U Imr g o Imr 



2^(56 + 8)1 (I) / V, 



J[2 



(27ri?o) 



fin)- 



{mi,m2) 

V, f{l/u) 



{2tiR^Y 



(3.28) 

2^(56 + 8), the 

number of massless boson/fermion pairs in the twisted sector, and e = — 1. The massless 



The contribution to the pressure is as in Eqs (3.15), (3.16) with n. 



''twist 
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bosonic content consists of 32 scalars in the (56,1) representation of Ei x SU{2) and 128 
scalars in the representation (1,2). The number of massless fermionic degrees of freedom 
follows by supersymmetry. 

Adding the contributions of the untwisted and twisted sectors together, we obtain 



Tin 



4 X 504 + 16 X 128 = 4064 



n 



V 



4 X 504 - 16 X 128 



-32. 



In addition, we have the choice with e = 1 in Eq. (3.16), giving 



n 



V 



(3.29) 
4064. We can 



also change Uy by considering Susy-breaking operators involving the right-moving gauge 



charges as in Eq. (2.48). 



3.4 Small mass scales from Wilson line deformations 

A generic supersymmetric heterotic background may contain in its spectrum massive super- 
multiplets whose mass is obtained by switching on non-trivial continuous Wilson-lines [27] 
- [29] . This is a stringy realization of the Higgs mechanism, breaking the initial gauge group 
G to a smaller one spontaneously. This statement is not absolutely correct for discrete 
Wilson lines corresponding to extended symmetry points where the gauge symmetry may 
enhance or even get modified. 

For our purposes, we restrict to arbitrary and small Wilson line deformations starting from 
a given supersymmetric background where Rj, I = 6, 7, . . . , 10 are of the order the string 
scale. This restriction ensures that the contributions to the thermal partition function of 
the momentum and winding states in these five internal directions will be exponentially 
suppressed in the limit where Rq and R^ are large. 

A systematic study of the effects of Wilson lines can be found in [27] - [29]. In the zero 
winding sector, a Wilson line just modifies the Kaluza-Klein momenta, and the corresponding 
Kaluza-Klein mass becomes 

{mi + y'} Q^f 

R] ^ R] ' ^^-^^^ 

where Qa is the charge operator associated to the Wilson-line y^. We distinguish three dif- 
ferent situations according to the direction J: 
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i) / = 5, where R5 is larger than the string scale. 

ii) / = 6, with i?e of the order of the string scale (with N = 4 and N = 2 initial supersym- 
metry) . 

iii) / = 6, . . . , 10, with Rj of the order of the string scale (with = 4 initial susy only). 

In the first case, / = 5, after a Poisson re-summation, the net modification to the partition 



function is obtained by the following replacement in Eq. (3.5): 



t2 



°° dt ^ 



— 2^ e "^0 



mi ,m2 



^2i7r{2m2+32) J/fQa 



t2 



]■ 



(3.31) 



mi ,m2 



The term in the brackets can be replaced with 



cos (27r(2m2 + 92) vlQa) = 1 -2sm^ (7r(2m2 + 92) vlQa) • 



(3.32) 



In the ii) and iii) cases, we can set the momentum and winding numbers to zero, mj = rij = 0, 
so that the extra modification in the partition function is the insertion of the term: 



(3.33) 



Incorporating the effects of the Wilson lines up to quadratic order, we get for the overall 
pressure: 



p ^ 2^4 ^ (2) 



+ M' 



£(§) 

3 

£(1) 

5 

712 



U 



(n^-2E.sin^(27rm2ygQ:) ) 
|(2mi + l)iu + 2m2p 



E 



mi,m2 



mi ,m2 



(M|-2E.M^in^(27rm2i/gg-)) 
|(2mi + l)iu + 2m2P 

- 2 Es sign(g) sin^ ((2m2 + l)7rylQi) 
\2miiu+ (2m2 + l)p 



r(i) ^ M^'^ - 2 sign(.)Mf sin^ ((2m2 + l)7rylQi) ^ ^^^^^^ 



3 



r7ii,m2 

In this expression, we have defined 



Mt 



\2miiu + (2m2 + 1)P 



47r^ i?? 



(3.35) 
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for the pair of boson/fermion states s and also introduced 

M| = 5^M2, = J]sign(s)M2, (3.36) 

s s 

where sign(s) indicates whether the state s contributes positively or negatively to riy and 
My , both being possibly negative. 

The following comments are in order: 

• In the above expression, y^Qa summarizes the contribution of the i?5-Wilson line in 
the term corresponding to the pair of boson/fermion states s. It does not introduce a new 
scale. 

• The Ms's introduce new mass scales in the theory, quahtatively different than T and 
M. The masses Mg are supersymmetric mass scales rather than susy-breaking scales like 
T and M. One must keep in mind that after minimization with respect to yj, physical 
scales like Mg are always proportional to an (infrared) renormalization group invariant scale 
like, for instance, Aqcd or the transmutation scale Qq, appearing in the electroweak gauge 
symmetry radiative breaking (of supersymmetric theories in the presence of soft breaking 
terms) [25] [26]. 

• The first two terms (which arise from the (91,92) — (1,0) sector) can be identified in 
the effective field theory as the thermal contribution to the pressure, Pthermai- Again the 
number is always positive being the number of the massless boson/fermion pairs in the 
initially supersymmetric background. This purely thermal piece is always positive. 

• The two last terms can be identified as minus the effective potential — Vcfr- This is 
naturally regularized in the infrared by the temperature scale T. This infrared regularization 
differs from that considered in [30], and used in [6], which is valid at zero temperature. 

• The number Uy can be either positive or negative depending on the way supersymmetry 
is broken. This shows that the sign of the one loop effective potential depends on the 
way supersymmetry is broken, as it can be seen in supergravity by utilizing super-trace 
arguments [31] [26]. 

3.5 Scaling properties of the thermal effective potential 

The final expression for P contains various mass scales: the two supersymmetry breaking 
scales which are the temperature T and the gravitino mass scale M, as well as the super- 
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symmetric masses Ms which are generated by the Wilson-hnes in the directions 6, 7, 8, 9, 10. 
The first identity follows immediately from the definition of P: 

which can be best seen by writing P as 

P = T^ p,{u) + p2{u), u^^, (3.38) 



where 

P4 



^ (F{u,yt) + F{u,yt)) , P2 = (^^^5") + G(u,yt)) , (3.39) 

-772 V / Tra V / 



71-2 V / 7r2 

and using the definitions 



^y^.y^) 2^ [(2mi + 1)2x^2 + 4mi]5/2' 



mi, 7712, s 



p/„ „a^ V" M^sign(g) cos(27r(2m2 + l)ygQa) . 
li^M^ cos(47rm2?/5(5a) 



mi ,m2 



^{^,y5) 2^ [(2mi + 1)2^/2 + 4mi]3/2' 

ti(s)M2cos(27r(2m2 
[4m2M2 + (2m2 + 1) 2] 3/2 



,,a^_ ^^sign(s)M2cos(27r(2m2 + l)y^,Ql) 



7711,7712,3 



Using standard thermodynamics identities, we can obtain the energy density p: 

rIP 

P = t"—-P^T' r^{u) + r2{u) (3.42) 

where 

r4 = 3p4 - up'^ , r2 = P2 - up2 (3.43) 

and the primes stand for derivatives with respect to u. In the sequel, we allow the Susy- 
breaking scales T and M to vary with time while fixing the super symmetric masses Mg, and 
investigate the back-reaction to the initially flat metric and moduli flelds. 
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4 Gravitational equations and critical solution 



We assume that the back-reacted space-time metric is homogeneous and isotropic, 

ds^ = -de + a{tf dnl H= (^^^, (4.1) 

where Qk denotes the three dimensional space with constant curvature k and H is the Hubble 
parameter. 

From the fact that —P plays the role of the effective potential and the relation between the 
gravitino mass scale M and the no-scale modulus $, 

we obtain the field equation for $: 

dP f dP\ 

^ + m^ = — = auy—j =-a{T\r^-?>p^)+T\r2-p2)). (4.2) 
We have made use of Eq. (3.43). 

For other fiat moduli y^j, with $ independent kinetic terms, the equation of motion is straight- 
forward to solve, 

y3, + 3i/0, = O =^ \^1 = %^ (4.3) 
where the Cj's are integration constants. 

Knowing the thermal effective potential —P, the energy density p as well as the field equation 
for the modulus $, we can derive the (one-loop) corrected space-time metric by solving the 
gravitational field equations. These are the Friedmann-Hubble equation. 



3H 



2 



2 2^"^' a 



and the equation that follows from varying with respect to the spatial components of the 
metric: 



2H + = - P - -^-^ - ■ (4.5) 



For our purposes, it will be useful to replace Eq. (4.5) by the linear sum of Eqs (4.4) and 



(4.5), so that the kinetic terms of $ and ipi drop out: 



9k 1 

H + 3H' = -- + -{p-P). (4.6) 



26 



4.1 Critical solution 



The fundamental ingredients in our analysis are the scaling properties of the thermal effective 
potential —P = — T^p4— T^p2- These scaling properties suggest to search for a solution where 
all varying mass scales of the system, M($), T and 1/a, remain proportional during time 
evolution: 

e"* = M(<l>) = — =^ H = -a^, M{<^)=uT, (4.7) 
7a 

with 7 and u fixed in time. Our aim is to determine the constants 7 and u. 



On the trajectory (4.7), the ^-equation is given by 



H + 3H'' = a'{ (r, - 3p,)— + (r, -P2)—] 



and the gravity Eq. (4.6) by 

1 1 

H + 3H^ = -2kj^M^ + (r^-p^) + -{r^ - ps)^ ■ (4.9) 

The compatibility of these two equations requires an identification of the coefficients of the 
monomials in M. The quartic terms give an equation for u, while the quadratic terms 
determine the sign of the parameter k and the magnitude of |/c7^|, 

= 2)^2 _ , ^^^4 = 4p4 , i.e. p4 + m = , for = , (4.10) 



- 2A:f = (r, - p,) \ , f-2A:f = i^^^ , for = ^ ) . (4.11) 



Eq. (4.10) reminds us of the equation of state for thermal radiation in five dimensions. In 
the absence of Wilson lines, where r2 = P2 = 0, we have that p = 4P, which is indeed the 
5d state equation for thermal radiation. 

When non trivial Wilson lines are turned on, their equations of motion fix them to be 
proportional to the renormalization group invariant scale Qq. Qq is supposed to be much 
smaller than M, Qo <^ M, so that the y^-terms can be consistently neglected. Within this 
approximation the relevant terms of the ?/j-effective potential are quadratic in y] modulo 
multiplicative terms arising from wave-function renormalization. If we define 



I, a 



Ri 
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we obtain, in a collectively qualitative way, the effective potential for each Wilson line (no 
sum over I and a): 



-Vj^ = Z^{yyRj,Qo)T' 



So we have that 



yl 
Ri 



d 



Z^{yyRi,Qo) = AUl-\og 



-V^ = 



Ri 



a,n2 



; 



(4.13) 



(4.14) 



diy-i/Ri) 

fixing the Wilson lines to be proportional to the transmutation scale Qq. We stress at this 
point that the are positive dimensionless constants, which, in general, are functions of the 
dimensionless ratio u = M/T. The u dependence of the is due to threshold corrections 
induced in the large scale regime: M,T ^ Qq. As we will show below, the u dependence of 
the plays a crucial role in our analysis. 



When all other equations of motion are satisfied, the gravity Eq. (4.6) is equivalent to the 
total energy-momentum conservation: 



d 
dt 



P 



l<i.2 

2 ' 2 







(4.15) 



This imposes an extra constraint when Wilson lines are switched on. On the critical trajec- 

(4.16) 



tory (4.7), and by using (4.14), one obtains : 



+ pf = 0, 



where we have defined the renormalized quantities 



E2? 



I,a 



yl 

Ri 



U 



P2 -UP2 



Ri 



(4.17) 



I,a 



to appear as well in Eqs (4.8), (4.9) and (4.11). 



Next, let us consider the Friedmann-Hubble Eq. (4.4) along the critical trajectory (4.7). It 
becomes 



6a' 



1 



6q;2 



3H' 



3k 



(4.18) 



The dilatation factor in front of 3H^ can be absorbed in the definition of k, Cr and Cm, once 



we take into account Eqs ( |4.20[ ) , ( |4.21[ ) and ( |4.22[ ) below: 

a* 
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(4.19) 



where, for cP' = 3/2, 



■^-2^ = ^P^ (4-20) 
^ ^4 = 7r:?-T4p4, (4.21) 



and 

Clearly, a necessary condition for the effective curvature k not to vanish is to have non trivial 



Wilson lines, at least in the direction 6, and to ensure Eq. (4.16). 



We note that Eq. (4.19) also controls the dynamics of a FRW universe, where space has 
constant curvature k and is formally filled with a thermal bath of radiation (since the sign of 
Cr can a priori be positive as well as negative). There can be an extra contribution, arising 
from the kinetic terms of some extra fiat moduli, that scales as 

Finally let us address a seeming puzzle. We said that along the critical trajectory and in 
the absence of Wilson lines, our thermodynamic quantities satisfy p = 4P. How can this 
situation correspond to a Ad universe filled with thermal radiation? The answer is that we 
must take into account the kinetic energy density of the modulus field $. When there are 
no Wilson lines and the kinetic terms of the other moduli are switched off, the Friedmann- 
Hubble equation gives $^ = (r4T^)/4 along the critical trajectory, and so the total energy 
density and pressure satisfy 

P.o.=VaT'+U'=^-^T\ (4.23) 
Thus ptot = 3Ptot which is the Ad equation of state for thermal radiation. 

4.2 Numerical study 

4.2.1 Without Wilson lines yl 

Let us consider the case where the Wilson lines along the Xs-direction are switched off. From 



Eq. (3.40), we obtain 

7r2 
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p,{u) = ^ fc/(n) + Kfiu)) , f{u) = u'f{l/u) , (4.24) 



thanks to the identities Tt^f{u) = F{u,yl = 0) and riyfiu) = F{u,y^ = 0). Since / + uf 



and / + uf vanish at the origin, Eq. (4.10) admits a universal solution u = for arbitrary 



Tiy in the range — < < n^. This solution corresponds to M{t) = at some finite 
T, and is associated to a 5-dimensional purely thermal system (i?5 = oo). The Ad effective 
description we have considered is not valid in this case. 



We are thus looking for non trivial solutions m > of Eq. (4.10 ), which we write in the form 

f + uf 



n 



V 



Tin 



(4.25) 



"T f + uf 

One can show that such non-trivial solutions (and consequently, non-trivial cosmological 
evolutions) exist only for models satisfying 



1 

- — < 



< 0, 



(4.26) 



The non vanishing root of Eq. (4.25) is an increasing function of the ratio ny/n^ satisfying 
u +00 when riy/n^ and u ^ when riy/n^ ^ —1/15, (See Fig. [l|. 




n 



V 



-0.06 -0.04 



■0.02 



Fi gure 1: The non trivial root u of Eq. U^Hm as a function of the ratio ny/rCp. 



The corresponding value of Cr in Eq. (4.21) is finite and always positive. This can be seen 
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by noting that the quantity n^f + riyf appearing in the expression for Cr equals —u{'rv^f' + 
riyf') > 0. The positivity follows since /' and /' are always negative and positive respectively. 

When Wilson lines in the Xg-direction are switched on, we can further study the sign of k. 
Defining G{u,y^ = 0) = M^g{u) and G{u,yl = 0) = m'^'^ g{u), the quantity in \A.n\ 
becomes 



3 

TX2 



MUu)g{u) + Ml'\u)~g{u 



g{u) = ug{l/u) 



(4.27) 



where M^(u) 



My\u) are obtained from M|. and My' under the replacement [yf/Ri) 



r(2) 



Ajc"('u)Qq. Since k and have same sign (see Eq. (4.20)), one finds that 



k>0 



u > 1 and — 1 < 



M^{u) 



< 



9[u) 



(4.28) 



(2) 

Note that a sufficient condition for k to be negative is that m < 1, or My > 0. For the 



particular value u = 1, Eq. (4.25) requires that riy/n^ ^ —0.0320. 

Let us detail the models we considered in the previous sections by computing the quantities 

(2) 



n 



y, n^, M|., My' to determine u and k. 



The r2 + P2 = constraint (4.16) ensures the existence of non trivial Wilson lines and 
imposes extra restrictions on the space of critical solutions. Supposing for simplicity that 
the Wilson lines are functions of u which are all proportional, Cj(m) 



critical solution with u determined by = Ap^ (Eq. (4.21): 
^ o d\og{g + sg) 







dlogu 



V 



where 



(2) 
V 



CM 



V 



one finds on the 



dlogc 
dlogu 



(4.29) 



s G [—1, 1] being a -u-independent number in this case. The above relation (4.29) reduces 
further the permitted values of s = My^'/M^. For rj small, one can show that (4.29) is 
satisfied when u is sufficiently large. The analysis for a generic value of rj is quite complex, 
and depends crucially on the particular choice of yj in the configuration space. In the 
explicit examples we will consider below, we will skip this analysis. However, this is not a 
real obstruction since, as we will see later, the simultaneous presence of non trivial Wilson 
lines in the directions 5 and 6 modifies the critical value of m in a highly non trivial way. 
The study of the //-constraint on the space of (n, ^5 ) goes beyond the scope of this work. 
Keeping that in mind we disregard the constraint when 1/5 = in the examples which are 
displayed below. 
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Model 1: 

We consider an = 4 heterotic string model, for example the i^g x theory on T^, at 
finite temperature and when supersymmetry is spontaneously broken by choosing Qj? = Qf- 
For this model 

n*^ = n*y = 2^x 504 . (4.30) 
Since riy > 0, there is no critical solution m > 0. 
Model 2: 

In the same Eg x E^ heterotic string model, the i?-symmetry operator used to break super- 



symmetry is chosen to he Qf + Qes + Q'es order to have < (see Eq. (2.43) for the 
T = case): 

n*T = 2^ X 504 ; n*y = 2^ x (-8) . (4.31) 

Then riy/n'^ = —1/63 ^ —0.0159, and so the model admits a critical solution u > 0. One 
finds numerically that u ~ 1.46 and 7^0^ — 441. 

Let us consider Wilson lines y'} in the direction 1 = 6 and define 



where a = 1,...,8 stand for the Cartan generators of the first Eg factor and a = 9, ...,16 
for the second. The derivation of the charges Qa of the initially massless states (see the 
Appendix) gives: 

M^ = 2'x 60 {{Y'f + ■ ■ ■ + (Y'y) , (4.33) 

and 

M^^ = 2^x (-4) {{Y^f + ■■■ + (Y^y) . (AM) 

Their ratio is Wilson line independent, My^ /M^ = —1/15, which is larger than —g{u)/g{u), 
and so the effective curvature is negative. 

Model 3: 



2 



In the N = 2 orbifold model (see section 3.3 for details), we set Qr = Qf + Qh + Qe7 with 



Qe^ being 1 for the spinorial representations of E^, decomposed in terms of 50(12) ones, 
and for the vectorial ones. We find the following: 

n*j, = 4064 ; n*y = 992 . (4.35) 
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Since nyjn^ is positive, this model doesn't admit a critical point m > at this stage. 
Model 4: 



This is the N = 2 orbifold model constructed in Sect. 3.3, with 



n*j, = 4064 , n*y = -32 . (4.36) 

Since ny/n"^ = —1/127 ^ —7.87 ■ 10~^, there is again a non trivial critical solution u > 0. 
Numerically, we find m ~ 1.90 and 7^0^ — 130. 

Before switching on Wilson lines in the direction / = 6, the gauge group is E7 x SU (2) x Eg. 
Let us consider arbitrary Wilson lines, Y"', a = 1, 16: Y^,...,Y^ for the Cartan generators 
of the 50(12) subalgebra of Ej, Y^ for the Cartan generator of the SU{2) subalgebra of E^, 
Y^ for the SU{2) factor and Y^,...,Y^^ for the Es- The mass square scales, computed in the 
Appendix, are: 

6 16 

M| = 624 ^{Y^y + 1248 (r^)^ + 736 {Y^f + 240 ^(F")^ (4.37) 

a=l a=9 

and 

6 16 

M^^ = -144 ^(F")^ - 288 (Y^ + 224 (Y^f + 240 ^(F'^)^ . (4.38) 

a=l a=9 

The smallest value of My'^ /M^ is reached when any F", a = 1,...,7, is non trivial (while 
= 0, a = 8, 16). In this case My^ /M^ = —3/13, which is larger than —g{u)/g{u), so 
that k is negative. 

4.2.2 With Wilson Lines yl 

Any model originally characterized by the quantities UT.rfy, M| and M^'^ can be deformed 



by switching on the Wilson lines 1/5,0 = 1,..., 16. We are looking for solutions of Eq. (4.10) 



written in terms of the functions defined in (3.40): 



F(u, yl) + uF\u, yl) + F(n, yl) + uF\u, yl) = . (4.39) 

We observe that the thermal contribution described by F{u,y'^) + uF'{u,y^) vanishes at 
u = and is bounded for large u. In all of the following examples, it is also positive. 
On the contrary, the "effective potential" corrections F{u,y^) + uF'{u,y^) vanish at m = 
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and diverge to +00 or —00 at infinity. Thus, in presence of arbitrary Wilson lines yf, the 
universal solution u — remains and we are looking for non trivial ones u > 0. These can 
only arise when F{u, yf) + uF'{u, yf) takes negative values. 

Model 1: 

Among the Wilson lines 1/5, (a = 1,...,8 for the first factor and a = 9, ...,16 for the 
second), wc choose to switch on either 

i) yl (with yl = ■ ■ ■ = yf = Q) OY 

ii) • • • = yl (with = . . . = = 0). 

In these cases and at fixed Wilson lines, F{u, yl) + uF'{u, yl) increases from to +00, and 
so there is still no solution w > 0. 

Model 2: 

For the Wilson lines defined in case ii) above, one finds that the critical solution -u > 
present before deformation is slightly shifted. However, in case i), the root m > is sent to 
+00 when y\ approaches the numerical value ~ 0.227 from below. For y\ above this critical 
bound, there is no non trivial solution anymore. 

Model 3: 

As said before, this model does not admit a critical solution w > when all Wilson fines in 
the fifth dimension are switched off. We consider four patterns of deformations, by switching 
on a Wilson line for a single Cartan generator: 

i) y\ (Cartan generator of 5*0(12) C E^) 

ii) yl (Cartan generator of SU (2) C Ej) 
Hi) (Cartan generator of SU{2)) 

iv) yf (Cartan generator of Es) 

In cases i) and iv), a large enough Wilson line deformation generates a non-trivial solution. 
We find the critical bound y^ ~ 0.895 for case i), and yf ~ 0.772 for case iv). In all these 
cases, the phase transition occurs when the limit at -u ^ 00 of the "effective potential" 
contribution F[u,yl) + uF'{u,yl) is switched from -\-oo to —00. In some sense, a solution 
u — +00 appears at the transition, and then decreases for a larger Wilson line, and finally 
reaches a non-zero minimal value. In cases ii) and Hi), since F{u, yl) +uF'{u, yl) is positive 
for any value of the Wilson line, the deformation does not create a non-trivial solution. 

Model 4: 
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We consider the same four patterns of Wilson lines. This time, the critical solution u > 
exists before deforming the model. We find that in the above cases, switching on a large 
enough Wilson line makes this solution disappear. However, two distinct behaviors are 
found: 

In cases i) and ii), we find again that the phase transition corresponds to a change in the 
behavior of the "effective potential" contribution at oo. In a similar way as before, the 
solution M > existing before deformation will increase when switching on the Wilson line, 
then go to +00 at the transition point, and disappear for a larger Wilson line. We find the 
critical bound yl ~ 0.105 for case i), and yl ~ 0.074 for case ii). 

For Hi) and iv), switching on the Wilson line makes the critical solution u > decrease 
towards zero, so that we are left only with the universal solution when the Wilson line is 
above some bound. In case Hi), this maximal value is yf ~ 0.205, while it is yf ~ 0.207 for 
case iv). 

Some remarks are in order 

• In all cases presented above, when a critical solution u > exists, the effective radiation 



term Cr given by Eq. (4.21) is strictly positive. 

• In some models, incorporating the Wilson lines in the fifth direction allows the critical 
solution M > to be close to or large. In other words, a hierarchy between the scales M 
and T can be found by tuning the moduli y^. 



5 Cosmological evolutions 

When a non degenerate solution u > exists, the model admits a well defined low energy 
description in four dimensions. The dynamics in this regime is controlled by the Friedmann- 
Hubble equation, whose behavior depends drastically on the signs of Cr and k. In all models 
we considered here, with initial = 4 or = 2 supersymmetry, Cr turns out to be positive. 
However, this situation may not be generic in more complex stringy examples with initial 
= 1 supersymmetry. For completeness, we briefly describe all possible cosmologies arising 
for any positive or negative value of c^. In addition, we allow non trivial time dependent 
profiles for the moduli 0j, giving Cm > 0. 
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5.1 Solutions for > 



• For k = 0, Cr > 0, Cm > 

When Wilson lines in the directions 6, 7, . . . , 10 are absent, the curvature k vanishes. In real 



time, the Friedmann-Hubble Eq. (4.19), 



3i7^ = ^ + ^, (5.1) 



can be used to express the time variable t as an integral function of the scale factor a, 

= a/^ r ^£2^ , (a>0), ao = ./^. (5.2) 

V Cr Jo + V Cr 

In this form, it is straightforward to draw a(t) (see Fig. |2]). The explicit solution is 

t{d) = if^o I T^-i/l + 72 ~ arcsinh ( -^J 1 , = ~JJ2 ' (^'^^ 




describing an expanding universe that starts with a big bang. As can be seen from Eq. (5.1 ), 

the slope d is infinite when a vanishes. At large t, the scale factor behaves as in the = 
particular case: 

a(t) = (l)'^'v^, (t>0). (5.4) 



Since the transformation t — — t is a symmetry of the Friedmann-Hubble equation, the 
previous expanding solutions have contracting counterparts and thus ending at t = with a 



big crunch. Finally, since the RHS of Eq. (5.1) is positive, there is no solution in Euclidean 
time. 

The effective field theory description always breaks down before the occurrence of a space-like 
singularity, when the temperature (in string frame) is of order the Hagedorn temperature. 
At this temperature scale, new stringy dynamics must be taken into consideration which can 
result into a phase transition, realizing the scenario of [16] . 

• For A; > 0, > 0, c„i > 

When we switch on Wilson lines we could have > 0. Rewriting the Friedmann-Hubble Eq. 



(4.19) in the form 



a± 



6k 



(5.5) 
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Figure 2: Cosmological evolution for the case k = 0, > 0, 



> 0. 



one expects a cosmological evolution satisfying < a < a+ should exist, while a solution 
with a scale factor greater than a+ should make sense in imaginary time only. In real time, 
one can actually express t as a function of a as follows 

(1 v^dv \ 

where 

1 r+ v^dv 

ti = ^ / , . 5.7 

From these expressions, one can see that the cosmological evolution starts with a big bang 
at t = tj. It expands until t = where the maximum size of the universe a_|_ is reached, and 
then it contracts until a big crunch occurs at t = — t,. 

To find a Euclidean solution, one needs to consider a scale factor greater than a^. It is 
then possible to find it by proceeding as before, or use the fact that such a solution can be 



obtained by analytic continuation of the expression (5.6) at t = 0, where a = a+. One can 
write 




_ , v^dv 1 v'^dv 

t = ±\ ti + — / —======== + 



kJo ^/{v^ + al){al-v'^) ^/lJa+ ^/-{v^ + al){v^ - al 



(5.8) 
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from which we derive 
T(aE 



as 



v^dv 



kJa+ v^(t;2 + a2)(i;2_a^) 
Fig. [3] represents the solutions a{t) and a£;(r). 



(5.9) 




Fi gure 3: Cosmological evolution for the case k > Q, Cr > Q, Cm > Q (in hold line). A big hang and a hig 
crunch are occurring at t = ti and t = —ti respectively. This solution is connected to a Euclidean one at 
t = —ir = that is asymptotically linear. 



In the particular case where Cm = 0, the solutions (5.6) and (5.9) are taking the explicit 
forms 



a(t) =a+^/l- ( ^ ) , iU<t<-U 




1 M 



fcV 3 



and 



whose shapes are similar to the generic case with > 0. 
• For k < 0, Cr > 0, Cm > 



This case is easier to deal with. Eq. (4.19) can be rewritten as 



[a^af' = |A:|(a^ + a^)(a^ + a^) , a± 



Cr ± y c'^ — I2\k\c„ 
6\k\ 



(5.10) 



(5.11) 



(5.12) 
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and admits the expanding solution 



t(a) 



1 



v'^dv 



^Jo V{v' + ai){v^ + ai)- 
(See Fig. Ill) After a big bang, the scale factor is growing linearly in time. 



(5.13) 




Fi gure 4: Cosmological evolution for the case fc < 0, > 0, > 0. 



The result for the particular case Cm = can be written more explicitly. The solution takes 
the form 



a{t) = 



t + to 



1 , (t > 0), 



; ^0 



to J — - " ' Y 3|^| ' - Y 3 

These cosmological solutions do not admit a sensible Euclidean continuation. 



(5.14) 



5.2 Exotic cosmologies with < 

Although in all explicit models we presented before Cr is positive, it is interesting to analyze 
the exotic situation with negative c^, which is not a priori forbidden in more general cases 
with = 1 initial supersymmetry. 

• For k = 0, Cr < 0, Cm > 

A cosmological evolution in real time only exists if is switched on. The scale factor 
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satisfies < a < ao where 




(5.15) 



Between a big bang at < and a big crunch at —ti > 0, a reaches a maximum oq at t = 0. 
At this time, an analytic continuation is allowed: A Euclidean solution satisfies > ao and 
goes to infinity for large positive or negative Euclidean time (see Fig. [s]). 




Fi gure 5: Cosmological evolution for the case k — Q, Cj. < Q, c,n > (in bold line). A big bang and a big 
crunch are occurring at t = ti and t = —ti respectively. This solution is connected to a Euclidean one at 
t = -iT = 0. 



• For /c > 0, < 0, Cm > 

Most of the considerations of this case are identical to the one derived for A; > 0, c^, > 0, 



> 0. In particular, the Friedmann- Hubble equation is still given by Eq. (5.5) and both 



the cosmological solution (5.6) and the Euclidean one (5.9) are valid. They are shown in 
Fig. [6] As long as > 0, the only qualitative difference with the case c.^ > is that 
the Euclidean solution has two symmetric infiexion points. However, when vanishes, the 
evolution in real time ceases to exist. 

• For k <{), Cr Cm > 

This case presents the most interesting features and involves either a first or second order 
phase transition in the early universe. The former case is the only one considered in this 
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Figure 6: Cosmological evolution for the case k > 0, < 0, c™ > (in bold line). A big bang and a big 
crunch are occurring at t = ti and t = —ti respectively. This solution is connected to a Euclidean one at 
t = —ir = that is asymptotically linear and has two symmetric inflexion points. 



paper, where a Euclidean solution has a finite action and thus can be interpreted as an 
instanton involved in a tunneling effect. These behaviors are qualitatively similar to the 
inflationary case studied in [6,12]. To be more specific we consider the Friedmann- Hubble 
equation in the form 

3(a^d) = 3|fc|a^ — |cr|a^ + Cm , (5.16) 
and discuss various regimes, depending on the value of the discriminant of the RHS 

5 = cl - I2\k\c^ . (5.17) 



i) When 5 > 0, there are two critical values for the scale factor 



a± 



\ 



Cr\±\lcl-I2\k\c„ 

Q\k\ 



(5.18) 



Eq. (5.16) then admits two distinct cosmological evolutions. The first one satisfies < a < 



a_ and corresponds to the usual dynamics between a big bang at tj < and a big crunch 
at —ti. The scale factor reaches a maximum a_ at t = 0. The second one describes an 
asymptotically linear contracting solution followed by an asymptotically linear expanding 
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one. The two branches are smoothly connected at t = 0, where the scale factor reaches a 
minimum value a+. Therefore, this solution is non-singular. (See Fig. ItI) 




Figure 7: There are two cosmological evolutions (in bold lines) for the case k < 0, Cj, < 0, > 0, when 
S > 0. The first one starts with a big bang at t ^ ti and ends with a big crunch at t ~ —ti. The second one 
has a contracting phase followed by an expanding one. These two branches are connected to each other by a 
first order phase transition via an instanton. 



The two cosmological evolutions are also related to one another by a double analytic con- 
tinuation: t = —IT and then t = Tf + it. Between r = and r = r^, a Euclidean solution 
whose action can be shown to be finite is allowed. It is thus an instanton between the two 
branches in real time and contributes to a first order phase transition. We note that when 
Cm = 0, the big bang / big crunch solution disappears since a_ vanishes. 



a) Let us turn now to the second case where the discriminant (5.17) satisfies 5 < 0. Eq. 



(5.16) does not admit any critical point and the scale factor is never stationary. There is a 
single cosmological evolution (and no Euclidean solution). It increases from a big bang at 
t = 0, while for large t, its time dependence becomes linear. Thus, close to the big bang, the 
evolution is similar to the first solution occurring when 5 > 0, while for large t its behavior 
is similar to the second expanding solution. Since there is an inflexion point at t = ti,-,f when 
a = \/2cml\cr\, the cosmological evolution for 5 < can be interpreted as a second order 
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phase transition between the same initial and final states encountered in the first order phase 
transition for 5 > 0. (See Fig. |8}) 




Fi gure 8: Cosmological evolution for the case A; < 0, < 0, Cm > 0, when (5 < 0. It describes a second 
order phase transition occurring at ti^j, between a phase that starts with a big bang to another phase that 
expands linearly in time. 



Hi) In the critical case 5 = 0, Eq. (5.16) admits two expanding cosmological evolutions 



which are asymptotic to a static one, a = Oq, where 

ao = a± = 4 / = flinf = Wfn^ , (5.19) 





together with two contracting ones obtained by time reversal. The first expanding solution 
starts with a big bang, while the second one is linear for large positive time. (See Fig. [91) 



6 Conclusions 

We have obtained several cosmological solutions in a large class of four dimensional heterotic 
string compactifications with spontaneously broken = 4 or = 2 space-time supersymme- 
try. The cosmological evolution is induced once radiative quantum and thermal corrections 
are taken into consideration. These corrections are calculated at the perturbative string 
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t 

Figure 9: There are two expanding (contracting) cosmological evolutions for the case k <Q, Cr Cm> 0, 
when (5 = 0. All are asymptotic to the static solution a = qq. 

level and shown to possess universal scaling properties. The reason is an underlying duality 
between the temperature and the supersymmetry breaking scale. 

Our solutions correspond to homogeneous and isotropic Friedmann-Robertson- Walker uni- 
verses. They are characterized by the ratio of the supersymmetry breaking scale to the tem- 
perature, and this ratio remains constant during time evolution. Even though Kaluza-Klein 
states associated to the supersymmetry breaking cycle are thermally excited, the equation 
of state governing cosmological evolution is identical to that of massless thermal radiation 
in four dimensions. This is due to the special relation between the no-scale modulus field 
associated to the supersymmetry breaking scale and the Hubble parameter: = 2H^/3. 
Universes with spherical, toroidal or hyperbolic spatial sections can be found once we incor- 
porate Wilson line deformations. 

In this paper we focused on the low temperature phase of the models. When the temperature 
is close to the Hagedorn temperature our effective field theory analysis breaks down and new 
stringy dynamics must be taken into consideration. It would be interesting to investigate 
if phase transitions can occur in these models as the temperature approaches the Hagedorn 
temperature, and whether such phase transitions result in non-singular time-dependent ge- 
ometries. To this extent it could prove useful to incorporate in our work the proposal of [23], 
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where such a phase transition is shown to occur in N — 4 heterotic string models, and study 
the cosmological imphcations. 

It would be interesting to extend our analysis to the N — 1 heterotic orbifold models, 
and for the cases where supersymmetry is broken spontaneously. In this class of models, one 
expects to find infiationary phases, once radiative and thermal corrections are properly taken 
into account. The analysis of [6, 12] reveals interesting transitions between such infiationary 
phases and radiation dominated phases with similar properties to those found in this work. In 
our examples, the coefficient of the term in the Friedmann-Hubble equation is positive. 
Perhaps among the = 1 examples it is possible to find models characterized by negative 
values of this coefficient. Then non trivial cosmological phenomena would occur, including 
first or second order phase transitions that allow for the possibility to realize the proposal 
for the creation of a universe from "nothing" [11] in string theory [6, 12]. 

The relation between the supersymmetry breaking scale with the temperature is a key prop- 
erty of our solutions. Suppose that such a scaling property persisted in an early universe 
epoch, and that initially supersymmetry was broken around the string scale. During such 
epoch, the Susy-breaking scale gets lower and lower as the universe expands and cools. At 
lower temperatures new dynamics may become relevant that can stabilize this scale. Such 
a scenario can give us a new perspective on how to handle the hierarchy and naturalness 
problems. 

Note added: A follow up of the present work can be found in [33]. The radiation era we 
have described is not only obtained for specific initial boundary conditions that select the 
critical trajectory. Instead, the radiation era is an attractor of the dynamics. It is reached 
asymptotically for generic initial boundary conditions. 

Acknowledgements 

We are grateful to Constantin Bachas, Ramy Brustein, Dieter Liist, Marios Petropoulos, Jan 
Troost and Fabio Zwirncr for useful discussions. N.T. thanks the Ecole Normale Superieure 
and CK. and H.P. the University of Cyprus for hospitality. 

The work of C.K. and H.P. is partially supported by the EU contract MRTN-CT-2004-005104 



45 



and the ANR (CNRS-USAR) contract 05-BLAN-0079-01 (01/12/05). N.T. and C.K. are 
supported by the EU contract MRTN-CT-2004-512194. H.R is also supported by the EU 
contracts MRTN-CT-2004-503369 and MEXT-CT-2003-509661, INTAS grant 03-51-6346, 
and CNRS PICS 2530, 3059 and 3747, while N.T. is also supported by an INTERREG IIIA 
Crete/Cyprus program. 



46 



Appendix A 



In this Appendix, we provide details on the derivation of the charges associated to the 
Cartan generators of the gauge groups of the heterotic string models we considered. We 



also derive the contributions to M|, and My-* (see Eq. (3.34)) associated to Wilson lines 
in the internal direction 6 only. Since the result is linear in the (2/g)^'s, the result for 
arbitrary Wilson lines in the directions 6, 7, 8, 9, 10 is obtained under the replacement 
[vlfK^^RD - E/(l/})V(4vri?|). 

Consider the = 4 heterotic model with gauge group £"8 x E^t,. When built out from 32 
worldsheet fermions as in the standard procedure, the vertex operators for the gauge fields 
associated to the Cartan generators of the first factor are given by 

^2i-1^2i^ ^ = 1,...,8, (A.l) 

whereas those associated to the Cartan generators of the second E^ are given by 

A2*-U'\ i = 9, . . . , 16. (A.2) 

In terms of representations of S'0(16), the adjoint representation of E^ decomposes as 

248 = 120 128, (A.3) 

where the 120 is the adjoint representation of 50(16) and the 128 is the spinorial represen- 
tation with positive chirality. The corresponding vertex operators can be written explicitly 
by bosonizing the 32 fermions into 16 free bosons iJj, i = 1, . . . , 16. This formalism has the 
advantage to make the roots of the Lie algebra appear in a clear way. For the 120 we obtain 

120: (eW±^^fc), {j^k)) ® idHj, j, G {1, . . . , 8}. (A.4) 

The 8 latter vertex operators correspond to the Cartan generators. For the 128 we have 

128 : e^''''^''^^"''''^^^"'""'"'"^*^**^ (A. 5) 

with the GSO constraint HLi = 1- go^-l here is to compute the quantity 

= ^ E (^al/6)^ (A.6) 

6 se248 
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where we denote by Q1 the charges of a state s with respect to the Cartan generators. 
We are going to consider the case of one gauge group; the generahzation to x E^ is 
straightforward . 

To begin with, we recall that the Cartan states are neutral. Then we are interested in the 
112 remaining states of the 120 adjoint representation. It is not hard to see that 

sel20 i^je{l,...,8} ei,e2=±l 

Therefore we get a quadratic polynomial in the |/g's. Noticing that this polynomial is invari- 
ant under the transformations — |/g and ^ yi, we obtain 

SG120 i=l 

Computing the {yiy term, we get a = 28. 
For the 128, we see that 

T^iQaVer- E \Uyl + ^2yl + --- + ifledyl) . (A.9) 

sel28 ei,...,e7=±l \ i=l / 

If we set yl = 0, the symmetries yl <-> —y^ and yl ^ y^, valid for i,j — 1, . . . ,7, guarantee 
that this polynomial will be of the form 

pi^ivir- (A.io) 

i=l 

Restoring y^^O gives a (|/g)^ term and crossed terms y^yg. However, y^ has been artificially 
isolated in the treatment of the GSO constraint: by isolating other 7/g's and using the same 
arguments, we can show that our polynomial is of the from 

E(^«^6f = /5E(^6r- 

sel28 i=l 

We obtain (3 = 32. 

It is then straightforward to evaluate the sums encountered before. We obtain 
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and when coupling the Scherk-Schwarz cycle to the helicity of the £^8 representation 



1=1 



In the N — 2 models, the orbifolding breaks ^ E-j x SU{2), under which the adjoint 
representation decomposes as 

248 ^ (133,1) ® (56,2) © (1,3). (A.14) 

The Cartan generators of Eg give the Cartan generators of E^ x SU{2): 

{idHi, idHe, i{dHy - dHs)) ; i{dHr + dH^). (A.15) 

Switching on arbitrary yl, . . . ,yl,yQ, we compute the charges of the various states step by 
step. 

In the 133, we have 7 neutral Cartan operators, 60 ladder operators in the Adj{S0{12)) 
sub algebra 

e'(^^^^H^\ e {1, . . . , 6}, (A. 16) 

2 ladders in the Adj{SU{2)) 

^±i{Hr-Hs)^ (A.17) 

and 64 ladders in a spinorial representation 

obeying a GSO condition. We see that yl has a particular role here. The latter states have 
charges ±| under the first six Cartan generators, and charges ±1 under the seventh. Using 
the same arguments as before, we see that the sum for the spinorial states is of the form 

(i^iyir + C^yiry (A.19) 

The polynomial we are looking for is therefore 

E {QM'-^oj^iyir + '^i'^yir + iG fE(^6r + (2y6')') =36^(y^)^ + 72(ye^)^ 

sel33 i=l \i=l / i=l 

(A.20) 
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Note that if we want to couple the Scherk-Schwarz cycle to the helicity of the Et, we have to 
compute also sign(s)((5ay6)^' where states in the spinorial representations of the S0{12) 
subgroup contribute with a minus sign. To get this sum we have to put a minus sign in front 
of the 64 part, so that 

E sign(.)(gX)' = 20^(y^)^ + 2(2^,)' - 16 [^(y^)^ + {2ylf 

sel33 i=l \i=l 

-^Y.{ylf-^Kylf- (A.21) 



For the (56,2) representation, we begin with the states with vertex operators 

The corresponding Qly^ are respectively 

±yi± {vl + yl) , ±yi± {yl - yl) . (A.23) 
Therefore the sum for these states equals 

T^^QaVer - 4 f^iy^r + 24 ivl + yir + 4 ^(^e)' + 24 {yl - ylf 

s i=l 1=1 

-8j](y^)^ + 48((|/eT + (l/6n- (A.24) 

i=l 

The remaining states to be considered have vertex operators 

gf(±i7l±-±H6±(i?7+H8))_ (A.25) 

For these, we get 



16(^E(^6)' + (2y6Tj- (A.26) 
Adding everything, we get the final result for the representation: 

E (^«^6)' = 24 j^ivD" + 48 (yir + 112 {yir. (A.27) 

se(56,2) 1=1 
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If we couple to the helicity, we also have 



E ^^Sn{smy^,r = -8 E(l/6)' + 48 {ylf - 16 {yl 

se(56,2) 



(A.28) 



i=l 



For the 3 of SU{2) the two states 



,±^(i^7+^^8) 



have charges ±2. So we get 



se(i,3) 



(A.29) 
(A.30) 



In the twisted sector, we encounter the representation (56,1), whose sum is obtained by 
switching off yf in the result obtained for the (56,2) representation and by dividing the 
result by 2. One then obtains, 

E = 12 E(l/6)' + 24 {yir, (A.31) 



se(56,l) 



1=1 



E sign(.)(gX)2 = -4 E(l/6)' + 24 (1/6^)1 (A.32) 

s6(56,l) i=l 

We also encounter the (1,2) representation, where the sum equals 2(?/g)^. 
Application to Models 3 and 4 

For model 4, we set Qr — Qf + Qh- If we consider Wilson lines corresponding to the 16 
Cartan generators of E^x E-j x SU{2), the result is 

16 6 6 

60 E(^6)' + 36 J2{yir + 72{ylf + 24 Y.{ylf + 48 {ylf + 112 {ylf + %{ylf 



Mr 



2,(2) 



i=l 



i=l 



"47ri?i 



32 12 E(l/6)' + 24 (i/D' + 128 {2{ylf) 



So we get 



/ 16 6 \ 

= 240 E(^6)' + 624 E(z/6)' + 1248 {ylf + 736 {ylf 

6 \ i=9 i=l / 



(2) 



47ri?; 



/ 16 6 \ 

77^ 240 Y.{ylf - 144 E(y^)' - 288 {ylf + 224 (y«)^ . (A.33) 



j=9 



=1 
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For model 3, we set Qr = Qa + Qh + Qet We get the same expression for M|,, while 



(2) 



16 



j=l 



1=1 



i=9 



32 -4 + 24 + 128 {2{ylf) 



i=l 



/ 6 16 

112 J] (1/^)2 - 800 {yir - 288 (i/e^)^ + 240 ^^(y, 



6)^ 



(A.34) 



i=l 



i=9 
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